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Abstract 

We consider fractional differential equations of order a £ (0, 1) for functions of one 
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is given. An analog of the Frobenius method for systems with regular singularity is 
developed. For a model example of an equation with a kind of an irregular singularity, 
a series for a formal solution is shown to be convergent for t > (if a is an irrational 
number poorly approximated by rational ones) but divergent in the distribution sense. 
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1 INTRODUCTION 



Fractional differential equations are widely used for modeling anomalous relaxation and diffu- 
sion phenomena; see [31 [H] for further references. Meanwhile the mathematical theory of such 
equations is still in its initial stage. In particular, a systematic development of the analytic 
theory of fractional differential equations with variable coefficients was initiated only recently, 
in the paper by Kilbas, Rivero, Rodriguez-Germa, and Trujillo [10] (see also Section 7.5 in [TT]). 
For equations of order a G (0, 1), of the form 

(D£ + u) (t) = a{t)u{t), t>0, (1) 

where Dq + is the Riemann-Liouville fractional derivative, or 

(D (q) m) (t) = a(t)u(t), t>Q, (2) 

where IB>( a ) is the Caputo-Dzhrbashyan fractional derivative, their main results are as follows. If 
a(t) = A(t a ), and A(z) is a real function possessing an absolutely convergent Taylor expansion 
on an interval \z\ < 8, then the equation (1) possesses a solution of the form 

oo 

U (t) =t a ~ 1 J2 a nt an , 0<t<#, 
n=0 

while the equation (2) has a solution 

oo 

U (t) = J2 b nt an , 0<t<6. 
n=0 

In both cases the solutions are unique, if appropriate initial conditions are prescribed. 

Thus, for example, the property of a-analyticity of the coefficient a(t) (defined above) 
implies a similar property of a solution of the equation (2). In fact, we have u(t) = U(t a ) where 
U is holomorphic in a disk {z G C, \z\ < 6}. The coefficient a may be complex-valued as well. 

The above results open the way for developing a theory of a-analytic solutions of fractional 
differential equations in the spirit of classical analytic theory of ordinary differential equations. 
Here we give some results in this direction. 

If in (2) a(t) = A(t a ) where A is an entire function, then the above results from [TD] with 
9 = oo (stated there in a weaker form, only for real arguments of analytic functions, than 
actually proved) show that the solution of the Cauchy problem for the equation (2) is of the 
form u(t) = U(t a ) where U is an entire function. Following [TQl [11] we call such solutions 
a-entire. In particular, that is true, if A is a polynomial. A natural question is about the order 
of U (here we investigate this subject just for the equation (2) since its properties are closer 
to those of ordinary differential equations). It is known (see, for example, [2] or [7]) that every 
nontrivial solution of the equation vS k \z) = A(z)u(z), k G N, with a polynomial coefficient A, 
is an entire function of order 1 + deg(A) / k. In this paper we prove that the orders of the entire 
functions U corresponding to solutions of (2) do not exceed (1 + deg(A))/a. As a — > 1, this 
agrees with the above differential equation result. On the other hand, if deg(A) = 0, that is 
A(z) = A, A G C, then U(z) = E a (Xz), where E a is the Mittag-Leffler function whose order is 
1/a [31 E], which shows the exactness of our general estimate. 
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Next, we investigate systems of fractional equations with regular singularity, that is the 
equations 

t a (Do+ u ) (t) = A{t a )u{t) (3) 

and 

t a u) (t) = A(t a )u(t) (4) 

where A(z) is a holomorphic matrix-function. Under some assumptions, we prove that formal 
power series solutions of (3) and (4) converge near the origin and develop an analog of the 
classical Frobenius method of finding a solution. For scalar equations, the latter problem was 
considered in [TT|, HI] . 

Finally, in order to clarify characteristic features of fractional equations with irregular sin- 
gularity, we study a model example, the equation 

t 2a (D^'b) (t) = Xu{t), AgC, (5) 

where, as before, < a < 1. Assuming that a is irrational and satisfies a Diophantine condition 
(that is a is poorly approximated by rational numbers), we construct a kind of a formal solution 
of (5) convergent for t > 0. We prove that the series for the formal solution does not converge 
in the distribution sense, within a theory of distributions associated with the fractional calculus 
(see [T5l ITT] ). Thus, the formal solution u(t) cannot be interpreted as a distribution solution. 
It is interesting that u(t) is closely connected with a class of analytic functions with irregular 
behavior introduced by Hardy [5]. 



2 Preliminaries 

2.1. Fractional derivatives and integrals [3|, 117] . Let a E (0, 1) be a fixed number. 
The Riemann-Liouville fractional integral of order a of a function (p 6 L\(0, T) is defined as 

t 

o 

The Riemann-Liouville fractional derivative of order a is given by the expression 

t 

(°o»( f ) = F(T^)l/( f - T )^M^ °<^ T - 



that is (Dq + v?) (t) = (/q+ c V) {t), provided the fractional integral /o+ c V is an absolutely 
continuous function. If ip is defined on the whole half-axis (0, oo), then Iq + <p and Dq + ip are also 
defined on (0, oo). Below we will consider just this case. 

Note that the Riemann-Liouville derivative is defined for some functions with a singularity 
at the origin. For example, if <p(t) = t d , d > —1, then 

(°^) w - fSr^/- (6) 
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so that Dq + ip = 0, if <p{t) —t al . For ip(t) = t , d > — 1, we have also 



r(d + i) _ fd+a 



T(ci+l + a) 



(7) 



The Riemann-Liouville fractional differentiation and integration are inverse to each other 
in the following sense. If <p> G L\(0, T), then Dq + /q+V 9 = V 9 - The equality Ig + Dq + tp = ip holds 
under the stronger assumption that ip = Iq + iP with some ip G Li(0, T). The latter is equivalent 
to the conditions of absolute continuity of /q+'V on [0)^] an d the equality (l ± a (p) (0) = 0. 

Let a function <£> be continuous on [0, T] and possess the Riemann-Liouville fractional deriva- 
tive of order a. The function 



i)( Q ) 



1 



ri-« 



t 

l J(t-T)-°<p{T)dT-t-*<p{0) 



1 



t a r(i - a 



■v(o) 



is called t/ie Caputo-Dzhrbashyan, or regularized, fractional derivative. If ip is absolutely con- 
tinuous on [0,T], then 



(D (Q) ip) (t) 



r(i - a) 



(t - r)- a ^\r) dr. 



(9) 



In contrast to Dq + , is defined only on continuous functions and vanishes on constant 

functions. In most of physical applications, equations with B)^ are used, because a solution 
of an equation with the Riemann-Liouville derivative typically has a singularity at the origin 
t = 0, so that the initial state of a system to be described by the equation is not defined. 

oo 

Let v(z) = ^2 c nZ n be an entire function. Consider the Caputo-Dzhrbashyan derivative of 

71=0 

the function 

oo 

l p(t)=v(t a ) = Y,C n t an . 

n=0 

It follows from (6) and (8) that 

oo 

(D< Q > cp) (t) = c n (3(n)t a ^ = (® a v) (t a ) 



n=l 



where 



and the operator 



T(na + 1) 
Y{na + 1 — a) 



(V„V) (z) = J2 C nP( n )z 



n— 1 



n=l 



(10) 



is known as the Gelfond-Leontiev generalized differentiation operator (see OH]; in fact, D a 
is defined for wider classes of functions). 
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The operator, right inverse to D a , has the form 



i 



(see Sect. 22.3 in |17j). It will be convenient to make a change of variables setting z = Re l9 , 
r = Rt a . Then 



R r I -i a— 1 

■ l/a 



Pa/) (^) = / P - 1 

o L 



f(re ie )dr. (11) 



If /(*) = £ f k z\ then 

fc=0 

and it is easy to check that 

@ a J) a v)(z) = v(z)-v(0). (12) 

2.2. A class of distributions. Spaces of test functions and distributions behaving rea- 
sonably under the action of fractional integration operators were introduced by Rubin [15J 
(for a brief exposition see also [17] ; both in [IT] and [15] there are references regarding other 
approaches and earlier publications in this field). Proceeding from [15] . it is easy to come to 
a class of distributions where the Caputo-Dzhrbashyan derivative is defined in a natural 
way. 

Let <S(0, oo) be the Schwartz space of smooth functions on [0, oo) with rapid decay at infinity. 
Denote 

S + = {(/?€ 5(0, oo) : y? (0 (0) = 0, / = 0, 1, 2, . . .}, 
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tp G S+ : / ^{x)x l ~ a ~ k dx = 0, k = l,2,...>, 



o 



OO 

= { cp e S + : J v {x)x a - l ~ k dx = 0, k = 0, 1, 2, . 
o 

These spaces are interpreted as topological vector spaces with the topologies induced from 
<S(0, oo); see [J5] for various descriptions of these topologies including the description by semi- 
norms. 

Together with the Riemann-Liouville fractional integration operator Iff + , it is convenient to 
use the operator 



(/» (t) = -L j{r - f)-V(T) dr, t>0. 
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If (p,if) are sufficiently good functions, for example, if (p G L p (0, oo), ^ G L q (0, oo), p,q > 1, 

11, 

— I — = 2 — a, then 

P 9 

oo oo 

y <^(x) (/o+» (a;) da; = y V(z) (a:) da; 



(see Section 2.5.1 in [IT]). We will write this in the notation 

IfV) = (V, - 

In particular, if ip — it', where it G 5(0, oo), then 

( V ,B^u) = (u',I 1 r\). (13) 

It is known [15] that J^~ Q acts continuously from $^_ _ct onto <3>" _1 . Therefore the identity (13) 
can be used to define u as a distribution from ($^ _Q ) , if u G C 1 (0, oo), and it' has no more 
than a power-like growth near zero and infinity. This definition agrees with the classical one: 
if u is continuously differentiable at the origin too, then 

oo oo 

1 



(u',li' a (p) = — r / u\x)dx / (t- x)- a <p(t)dt 

N ' r(l — a) ' ' 



o x 

oo 



— - — - I <p(t) dt lit- xY a u'(x) dx = (if, D (Q) u) 
Til — a) J I 



o o 



where IJW u is understood in the sense of (9). 

A typical example of a function from is the function 



x a (x) = x exp I — I sin ^— logs ) . 

It is clear that K a G S + . Next, if 

( \ ( lo & 2x \ ■ f 71 , 
x(x) = exp I — I sin ^— log a; 

then we can write explicitly the Mellin transform 

oo 

k{z) — [ x z ~ x x(x)dx. 
Namely, by the formula (4.133.1) from [4], 

oo 







x(z) = 2 J e"' 2/4 sinh(zt) sin y dt = 2^e z2 ~^ sin nz. (14) 



In particular, 

oo 

J >c a {x)x l - a - k dx = x(-k) = 0, k = 0, 1, 2, . . . , 
o 

so that indeed x a e 

In order to have a full concept of a class of distributions, one needs a result regarding 
density of a space of test functions in some space of integrable functions. This gives a one- 
to-one correspondence between ordinary functions and distributions they generate. Here we 
present such a result though it will not be used directly in this paper. For similar properties in 
other situations see [To] . 

Proposition 1. The space is dense in L p ((0, oo), t~ l dt) , 1 < p < oo. 



Proof. Let R = J H a (i)t- l dt. We have 



Denote 



Then 



R = x(a - 2) = 2v^Fexp ^(q - 2) 2 - ^ sinvra > 0. 



**(a0 = ^a{x N ), N = 1, 2, ... . 



zn(x)x 1 dx 



R 



X, 



,{t)r l dt = i. 



(15) 



Suppose that / G L p ((0, oo),i Mt). Consider the so-called Mellin convolution 



{z N *m /) (t) = / z N {r)f - r-hlr 



that is actually the convolution on the multiplicative group (0, oo) (note that t 1 dt is a Haar 
measure on that group). Obviously, *m f £ ^ 1 i T a - Denote by || ■ \\ p the norm in L p ((0, oo), 
t^dt). 

Using (15) we can write 



(ZN *M /) (*) " /(*) = / Z N (T) 

By the generalized Minkowski inequality, 



/ 



/(*) 



"Mr. 



\z N *m f - f\L < / \z n (t 



f 



f(t) 



1/p 

r x dt \ dr. 
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Next we use the L p - continuity of shifts on the multiplicative group (0, oo) (see [9], Theorem 
20.4). For any e > 0, there exists 5 > 0, such that 



/(-)-/(*) 



i/p 

t~ L dt \ <s, 



if It — 11 < 5. Thus, 



\z N *Mf-f\\ p <2\\f\\ p J \z N (r)\T- L dr + e J \z N (r)\r- L dr. 

|t-1|>4 |t-1|<5 



Note that 



oo oo oo 

J \z N {T)\T~ x dr = ^ J \>c a {x N )\ x~ l dx = ^ j \x a (t)\r x dt = G x 



where the constant C\ does not depend on N. On the other hand, 

T2 1„„2 



\z N (x)\ < C 2 Nx N{a ~ 2) 



exp 



N 2 log z x 



If \x — 1| > 8, then log x > b > 0, so that 



(16) 



\Z N {X)\X 



-^dx<C 3 e-^(l + 8) N ^^0, 



1+5 



as A" — > oo, and (for 5 < 1) 

l-<5 log(l-<5) 



y \z N (x)\x~ l dx < C 2 N J exp jiVt(a - 2) 



N 2 t 2 



dt 



ATlog(l-<5) 



C 2 y exp |s(a - 2) - i-j -> 0, 



as A" — > oo. 

As a result, we see that, if A" is large enough, the first summand in (16) does not exceed 
2||/|| p £, while the second < C\S. Thus, 



\zN*Mf-f\L — >0, as Af^ oo. 



2.3. On ratios of the Gamma functions. We will often use the function 

r(t+i) 

*> = r( t +i- n) - "K^- < 17 > 
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Here we collect some of its properties. 
If t > a — 1, the integral representation 



oc 



1 f e --rt e («-D- (i- e -') a - l ds 



p(t) r(a) 
o 

holds (see Chapter 4 in |13J). It follows from (18) that the function t i— > — - r^r is strictly monotone 

decreasing and — -r^- — > 0, as t — > oo. Since r(t + 1 — a) has a pole at £ = a — 1, it is seen from 

(17) that p(t) -> 0, as t -> a - 1 + 0. 

If— 1 < i < a — 1, then (by a well-known identity for the Gamma function) 

T(-t + a) sin7r(t + l) 



Pit) 



T(-t) sin7r(t+l-«) 



The integral representation for the ratio of the Gamma functions [13J leads, after an elementary 
investigation, to the conclusion that p(t) is strictly monotone increasing from — oo to 0. 

Thus, we conclude that on the interval (—1, oo) the function p(t) is strictly monotone 
increasing from — oo to oo. The inverse function 7(A) solving the equation 

p{t) = A, A G R, 

is a well-defined continuous function. Note that p(0) = — -, so that 7(A) > 0, if A > 

r(l — a) 

1 



T(l -a)' 

It is known ([15]. Chapter 4) that 

^±1^^(1 + 0(0), t^oc, (19) 

if b > a. In particular, 

pit) ~ t a , as t — >• 00. 

For the sequence j3(n) defined in (10) and appearing in the definition of the Gelfond-Leontiev 
generalized differentiation operator D a , we have f3(n) = p(an), so that the above asymptotics 
implies the relation 

(3(n)~Cn a , 71 = 0,1,2,.... (20) 



3 a-Entire Solutions 

Let us consider a-entire solutions of the equation (2) with a(t) = A(t a ) where A is a polynomial 
of degree m > 0. We assume the initial condition u(0) = uo- 

Theorem 1. Under the above assumptions, the solution u(t) of the equation (2) has the form 
u{t) = v(t a ) where v is an entire function whose order does not exceed (1 + m)/a. 
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Proof. Seeking the function v, we have (T) a v) (z) = a(z)v(z). Let us apply the operator 3 C 
(see (11)) to both sides of this equality. We get from (11) and (12) that 



R 



v(Re iy ) - v(0) 



which implies the inequality 



aT(a) 



l/a 



- 1 



a-l 



a(re lti )v(re l0 ) dr. 



\v{Re ie )\ < |«(0)| + C 



l/a 



a-l 



\a(re ie )\\v{re l9 )\dr 



(here and below we denote by the same letter C various positive constants). 
We have the asymptotic relations 



t« - -(t-l), as t^l + 0; 
a 



Therefore 



so that 



a-l 



a-l 



a-l 

~ t <* , as t — > oo. 



t*-l) <C(t-l) a -H 



a-l 



, t>l, 



R 



\v(Re ty )\ < \v(0)\ + CR 2 -«~ a J {R-r) a '\^\a{re l0 )\\v{re lt, )\dr. 

o 

Since a(z) is a polynomial of degree m, we get 

R 

\v{Re ie )\ < |v(0)| + CR 2 -^~ a [ {R-r) a - l r«- 1+m \v{re i6 )\dr 



where C does not depend on R, 9. Fixing 6 and denoting 

\v{re %e )\ 



w(r) 



_2-i-Q 
f a 



we come to the inequality 



R 

w{R) < |t;(0)l + C [{R- r) a - 1 r 1 - a+m «;(r) dr. 
R 2 c a J 



(21) 



Now we are in a position to apply Henry's theorem (see Lemma 7.1.2 from [8j), which states 
that the inequality (21) implies the inequality 



W ( R "> - 2 _l_ Q £a,2-a+m(CR) 

R a a 
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where £ a ,a{s) is a certain function admitting the estimate 

£ a A s )< Cs ' A a J exp — ■ -s * 

\a + a — 1 

Thus, 

\v{Re i6 )\ <Cexp (/ii? 1 ^) 
for some fi > 0, as desired. ■ 

4 Regular Singularity 

4.1. Formal and a-analytic solutions. Let us consider systems of equations of the form 

t a (D£ + u) {t) = A(t a )u(t) (22) 

where 

oo 

A(z) = A + J2A m z m , 

m=l 

A m are n x n complex matrices and 

\\A m \\<Mfx m (fi>0), m = 0,1,2,.... 
Suppose we have a formal series 



u(t) = J2u k t a \ u k eC n . (23) 

k=0 

Let us substitute the series (23) formally into (22). We get, in accordance with (6), that 

oo oo 
k=0 m,k=0 

where (3(k) is the sequence (10). Collecting and comparing the terms we find that 

i 

(3{l) Ul = J2 A kUi- k , 1 = 0,1,2,..., 

k=0 

or, equivalently, 

A u = — -u ; (24) 

1 (1 — a) 

i 

[A -(3(l)] Ul = -J2 A kUi-k, l>l. (25) 

k=l 

It is natural to call the formal series (23) a formal solution of the system (22) if the relations 
(24), (25) hold. 
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Proposition 2. If a formal series (23) is a formal solution of the system (22), then the series 
(23) is absolutely convergent on some neighbourhood of the origin. 

Proof. It follows from (19) that 

H [v4 - /3(0] _1 II < l>l . 
In particular, we may assume that 

ii [A) - /^or 1 ii < i, i>io- 

Considering, if necessary, \u(t) instead of u(t), with |A| small enough, we may assume that 
IKH < 1- 

Let us choose so big r > that < r l for / < l and 



oo , 

fc=i 



Then 



\\ui\\ < r l for all I. 
Indeed, if this inequality is proved up to some value of I > l , then 

i+i 

Mr l+l 



\ u i+i\\ < 



i+i 

A k u l+ i_ k 



k=l 



l+l 



k^+l-k 



< r 



2+1 



k=l k=\ 

and the above inequality implying local convergence in (23) has been proved. 
The above arguments remain valid for systems of the form 

t a (B^u) (t) =A(t a )u(t). 



(26) 



The only difference is that, instead of (24), we get the relation A Q u = 0, just as in the classical 
case (see [B]). 



4.2. Model scalar equations. Consider the equation 



(27) 



By the relation (6), a solution of the equation (27) is const -t 7 ^ where 7 (the inverse function 
to p) was defined in Section 2.3. For example, if A = 0, then we have 7(0) = a — 1. 

If we consider an equation similar to (27), but with the Caputo-Dzhrbashyan derivative, 
that is 

t a (B^ip) (f) = \<p(t), (28) 

then the constant function is a solution of (28) for A = 0. Suppose that A 7^ 0, and if is a 
solution of (28), that is <p G C[0, T], the function 
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is absolutely continuous, and (28) is satisfied with 

*) (*) = j t M - 



We have 

1 

(JjfV) (t) = J (l- s)-Mst) da, 



and since <p is continuous, 

(/o+» (t) — ► 0, as t -> +0. 

It is known [T7] that in these circumstances Iq + Dq + (p — (p. Note also that 7q+ transforms 
the function t~ a into the constant T(l — a). Dividing the equation (28) by t a and applying Jq + 
to both sides, we find that 

t i 
<p{t) - <p(0) =*-f(t- r)- 1+a T->(r) rfr = A_ |(1 - s)" 1 ^- V^) 



Passing to the limit, as t — > 0, and taking into account the continuity of (p, we obtain the 
identity 

*jg /<!-.)-■".-.* = 0. 



whence y?(0) = 0. 

Thus, for A 7^ 0, the equation (28) is equivalent to (27), if (27) is considered for continuous 
functions vanishing at the origin. The power solution Ct" 1 ^ belongs to this class, if A > fnz^y- 
It may be instructive to see these solutions, satisfying the equation 

t 

^ ^ ■ - -l+a_-a, 



as examples of non-uniqueness of solutions of linear Volterra integral equations occurring due 
to the singularity of a kernel. 

4.3. Systems with good spectrum. Let us consider the equation (3) with A(z) = 

oo 

A m z m where A m are complex n x n matrices, the matrix Aq is Hermitian, and the series 

m=0 

converges on a neighbourhood of the origin. Without restricting generality, we may assume 
that 

A Q = diag(Ai,...,A„), Ai,...,A n eR. 

Following the classical method (see, for example, pQ) we look for a matrix- valued solution 
{a fundamental solution) of the equation (3), in the form 

u{t) = S{t a )tfj{t) (29) 
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where ip{t) = diag (t^ Xl \ . . . , t 7 ( An )), S(z) = Yl a v zV '> °v ( u — 1) are some unknown matrices, 

(T = I- 
We have 



= J2 <?» diag . . . ; t 7(A,0+a^ ^ 



Ml. , 

whence 

oo 

t a (D° + u) (t) = ^ diag (t^ +au , p^ +av ) (30) 

where 



R v u = diag 
On the other hand, 



r( 7 (Ai) + ai/ + 1) r( 7 (A n ) + az/ + l) 



r( 7 (Ai) + cm/ + 1 - a) ' " ' ' r( 7 (A„) + au + 1 - a] 



OO / V 



A(t a )u(t) = J2 [J2 A rn°»-m diag ( t ^+™ , . (31) 



u=0 \m=0 



Note that Ro = diag(Ai, . . . , A n ) = A , and since cr = I, the coefficients corresponding to 
v = in (30) and (31) coincide. Comparing the rest of the coefficients, we obtain the following 
system of equations for the matrices a k . 



k-l 



a k R k - A a k = ^ A k -iai, k > 1. (32) 



1=0 



For each k, the matrix equation (32) for a k has a unique solution if the spectra of the 
matrices R k and Aq are disjoint (see Appendix A.l in [1]), that is 

r( 7 (A,) +«* + !) 



r( 7 (Ai) + + 1 - a) 

for all i,j G {1, . . . , n}, or, equivalently, since the left-hand side of (33) equals p( 7 (Aj) + ak), 

7 (Aj) - 7 (Ai) 7^ afc, for all i,j G {1, . . . ,n}. 

We call our system (3) a system with good spectrum, if 

7 (A,) - 7 (A 4 ) i aN, for all i,j G {1, . . . , n}. (34) 

This definition extends the classical one [T], since for a = 1 we would have p(i) = 7 (i) = t, and 
the condition (34) would mean that the eigenvalues of A must not differ by a natural number. 

Theorem 2. If a system (3) has a good spectrum, then it possesses a fundamental solution 
(29) where the series for S(z) has a positive radius of convergence. 
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Proof. By the asymptotic relation (19), 



r( 7 (A i ) + ak + 1) 



(ak) a (l + 0(k~ 1 )) , k 



r(7(Aj) + ak + 1 - a) 
for all j = 1, . . . , n. Therefore 

(ak)~ a R k = I + 0(k- 1 ), k-^oo. 



oo, 



(35) 



Let us divide both sides of the equation (32) by (ak) a . The resulting equation, considered 
as a system of scalar equations for n 2 elements of the matrix o~ k , has the coefficients bounded 
in k and the determinant, which is different from zero for each k and tends to 1, as k — > oo. 
This implies the estimate 



\a k \\ < ak 



k-l 

A k -i<yi 



1=0 



k > 1, 



(36) 



where the constant a > does not depend on k. It follows from (35), (36), and the convergence 
near the origin of the power series for A(z) that 



fc-i 



\0~k\ 



< a x k- a b k ~ l \\ai\\ , k>l, 



1=0 



where a\ and b are positive constants independent of k. 

Define a sequence {sfc}^° of positive numbers, setting s = 1, 



fc-i 



s k = ai k- a Y, hk ' 1 ^ k > 1 - 



1=0 

The induction on k yields the inequality ||o"fc|| < s k for all k > 0. On the other hand, 

k 



s k+1 = a 1 (k + l)- a J2 bk+1 ~ l 



si 



1=0 



(k+l)- a 
k~ a 



k-l 



a i 



k' a b^b k - l Sl + bs k 



=0 



(fc + l)-° /, . , _ Q , x (k + l)~ a b , , 

(bs k + ai k bs k ) = — {l + ai k ) 



Therefore 



k~ a 



as k — > oo. 



Sk- 



(37) 



It follows from (37) that the series s k zk has the convergence radius b 1 (see Section 2.6 



k=0 



in [20]). Moreover, the series a v z u converges for \z\ < b 1 . 



u=0 



For the equation (26), a similar construction is valid, if we assume that Ai,...,A n > 
1 

T(l -a)' 
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5 Irregular Singularity: an Example 



5.1. A formal solution. In this section we construct a solution, in a sense to be specified, of 
the equation (5). Looking at classical first order equations, corresponding formally to a = 1, we 
have to consider the equation t 2 y'(t) = Xy(t) whose solution is y(t) = exp(— At -1 ). Therefore 
it is natural to seek a solution of the equation (5) in the form 

CO 

u(t) = Y,c n t~ na , c n eC. (38) 

n=0 

A fractional derivative of any term in (38) with n > a -1 does not make sense classically. 
However we may apply the distribution theory from Sect. 2.2. Below we understand the 
fractional derivative U>^ in the sense of (13). 

Proposition 3. (i) If pi < 0, ji ^ —1, —2, . . then 



r(/i + l-a) 



(39) 



(ii) If k is a natural number, then 



D (a) t' k = } 1 -r k - a log t. (40) 

(k - l)!r(-fc + 1 - a) B V ; 

Proof. Let up E $^T a . By (13), 

(» (a) f, <p(t)) = ti (f" 1 , (t)) . (41) 

It is clear that the right-hand side of (41) is an entire function of \x. For /i > 0, by virtue of (7), 

<f-\ (/!-» (t)> = (/fcT- 1 ,^)) = T{ ^_ a) (*™Mt)) ■ 

For /j, < 0, /j, 7^ —1, —2, . . ., the analytic continuation gives the equality (39). 
Next, consider the entire function 

Note that F(—k) = 0, k G N, by the definition of the space We have 

F'(n) = (logt-e^-^ loet ,uj). 

In particular, 

CO 

F'(-k) = J r a - k iogt-ip{t)dt. 
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As fj, belongs to a small neighbourhood of the point —k, = F'(—k)(fi + k) + o(/i + k). 

(_l)fc-i 

Since the residue of + 1) at fj, — —k equals — — (see Section 4.4.1 of [IS]), we see that 

(k- 1)! 

the function + l)F(n) is holomorphic at [i = —k (in fact, it is entire) and 



oo 

k-1 





Now the equality (39) implies (40). ■ 

Returning to (38), we will formally apply B)^ termwise and find the coefficients c n compar- 
ing resulting terms in (5). It is clear from the equality (40) that such a procedure would fail 
if some of the numbers na are integers. Thus, we have to assume that a is irrational. Using 
(39) we find that (formally) 

n=l v ' 

Substituting this into (5) we come to the recurrence relation 

r(l - (n + 2)a) 

c n+i = A— — — v ^c„, n > 0, 

T(l — (n + l)a) 

and it is easy to find by induction that 

rq-fo + i)^ 
Cn ~ x r(i-a) c °- 

Thus, we have found a formal solution 

oo 

u(t) = — ?L-r r(l - (n + l)«)t- na (42) 

1 ^ n =0 

of the equation (5). Using the identity 

r(z)rd-z) 



sin TTZ 

we can rewrite (42) in the form 

oo 

C 7T 



= r(l - a) ^ sin(vr(n + l)a) T((n + l)a) ' ( } 



5.2. The convergence problem. The convergence of the series (43) depends on the 
arithmetic properties of the irrational number a. It was shown by Hardy [5] that a can be 
chosen in such a way (to be well approximated by rational numbers) that the series (43) would 
diverge for small values of t. 
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An irrational number a G (0, 1) is said to be poorly approximated by rational numbers if 



p 

there exist such e > 0, c > that for any rational number -, p, q G N, 

9 



P 



> cq 



-2-e 



(44) 



By the Thue-Siegel-Roth theorem (see [18]) such are all algebraic numbers. 

The first statement of the next theorem is actually contained already in the paper [5]. 

Theorem 3. If a is poorly approximated by rational numbers, then the series (43) converges 
for any t > 0. However this series diverges in the space of distributions ($ 1 t r Q ) ■ 



Proof. It follows from (44) that 



\qa — p\ > cq 



-l-e 



for any p G N, 



so that 



dist((ra + l)a, Z+) > c{n + 1) 



-l-e 



and taking I G Z + , such that \{n + l)a — Z| < — , we find that 

|sin(7r(n + l)a)| = | sin(vr((n + l)a-/)| > 2\(n + l)a-l\ > 2 dist((n+ l)a, Z+) > 2c(n + l) 



-l-e 



Using the Stirling formula we obtain that the series (43) converges for each t > 0. 

To prove the second assertion, consider (t~ an , x a (t)) where the function x a was defined in 
Sect. 2.2. We have, by (14), that 



x a {t)) = J t~ an x a {t) dt = Kail — an) = Tk{a — 1 — an) 
o 

= 2y/ixex\) (^i a — 1 — otn) 2 — —^j sin(7r(a; — 1 — an)). 

Now we can give a lower estimate of the coefficients in the series (43) understood in the 
distribution sense: there are such a, a\, b > that 



1 



(r 



.(*)> 



> a(n + l)~ b(n+1) | sin(vra(n - 1))| exp(a 2 n 2 - 2a(a - l)n) 
sin(7r(n + l)a) T({n+l)a) ~ v ; I v v ))\ v\ \ ) j 

> ai {n + i)-Kn+i)-i-e eX p(a 2 n 2 - 2a[a - l)n) 

= a\ exp {a 2 n 2 — 2a(a — l)n — (6n + b + 1 + e) log(n + 1) } — > oo, 



as n ^ oo. Therefore the series (43) does not converge in a V. 
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